Dissociation of Feshbach Molecules into Different Partial Waves 
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I. INTRODUCTION 

The association of ultracold molecules from ultracold 
atomic gases using Feshbach resonances was a major re- 
cent breakthrough in the field of cold molecules [1-9] . Ex- 
perimentally, the method requires a slow magnetic-field 
ramp that crosses a Feshbach resonance in the proper 
direction. The molecules can be dissociated back into 
unbound atom pairs by ramping the magnetic field back 
through the Feshbach resonance. Initially, the atomic gas 
is so cold that only s-wave collisions are relevant in the 
gas. The molecular bound state that causes the Fesh- 
bach resonance is usually also an s-wave state, and the 
outgoing wave in the dissociation is again an s wave. 

If one chooses a Feshbach resonance where the molec- 
ular bound state is not an s-wave state, then one might 
wonder if outgoing waves other than the s wave can be 
produced. Of course, if the outgoing wave is as cold as 
the incoming one, the s wave will dominate again. But 
if the magnetic-field ramp for the dissociation is fast, 
then kinetic energy can be added during the dissocia- 
tion [10, 11]. It thus seems feasible to populate outgoing 
higher partial waves. This prompts many questions: If 
one creates a higher partial wave, will there still be an 
s-wave component? Do the different partial waves form a 
coherent superposition or an incoherent mixture? What 
determines the relative phase and the branching ratio? 
And how fast is the dissociation process? Are certain 
Feshbach resonances better suited than others for creat- 
ing a large fraction of a specific outgoing partial wave? 
Answering these questions is nontrivial and no theory has 
been developed on the subject, yet. 

The key to a theoretical description of the dissociation 
process lies in the observation that the dissociation is 
"half a collision". In a full collision, two atoms come 
together and then separate again. In the association 
and dissociation using Feshbach resonances, the exper- 



imenter can "freeze" the population in the middle, after 
the atoms came together. He can even choose, how much 
time is spent between association and dissociation. Still, 
the association and dissociation can be regarded as the 
first and second half of one collision. The concept of a 
half collision proved useful in other contexts before (see 
e.g. Refs. [12-15]). Unfortunately, the models developed 
there are not directly applicable here. The question is 
then how the dissociation can be linked quantitatively to 
scattering theory. The objective of the present paper is 
to establish this link. 

The motivation for this investigation comes from an ex- 
periment we performed recently, where dissociation into 
two partial waves, s and d, is observed [16]. The experi- 
ment employs *^Rb, where the dissociation is particularly 
interesting because of the presence of a d-wave shape res- 
onance. Our theoretical studies in the present paper are 
geared towards explaining the results of this experiment. 
A very brief summary of the model was already presented 
in Ref. [16]. 

The paper is outlined as follows: Section II begins with 
a brief summary of some basics of scattering theory. Sec- 
tion III presents the theory of scattering resonances for a 
single partial wave. The d-wave shape resonance in *^Rb 
is introduced as a specific example. In Sec. IV, mag- 
netically tunable Feshbach resonances in the collision of 
ultracold atoms are introduced. Section V describes some 
basics of scattering resonances with more than one par- 
tial wave and then discusses the combination of the shape 
resonance and the Feshbach resonance. The stage is then 
set for Sec. VI, where the link between scattering theory 
and the dissociation of molecules is presented. The decay 
rate, branching ratio, and the relative phase between the 
partial waves as observed in the experiment of Ref. [16] 
are explained with this model. 
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II. BASICS OF SCATTERING THEORY 
A. General 

The problem of scattering two particles off one another 
is easily separated into center-of-mass and relative coor- 
dinates. The ccntor-of-mass motion is trivial so that the 
problem is in the relative motion. The latter is equiv- 
alent to the scattering of one particle with the reduced 
mass rrired off a potential V{r). The relative motion is 
characterized by the wave vector k. The corresponding 
kinetic energy is E = T^k'^ /{2mred)- 

Scattering theory is usually formulated as a time- 
independent process with an incoming plane wave with 
k pointing along the i;-direction. One can show that 
the scattered wave at large radius falls off radially like 
a spherical wave. Hence, the scattering wave function 
has the asymptotic form 

^(+)(,-^ — e^fc^-^mc^)— , (1) 

r 

where r,'d,ip are spherical coordinates. The scattering 
amplitude /(i?, (p) is related to the differential scattering 
cross section 

^ = l/Wrf, (2) 

where dO, = sm-d M dLp is the differential solid angle. Fi- 
nally, the total cross section a is obtained by integration 
of the differential cross section over the full solid angle. 
The task in scattering theory is to determine /(^, ^p) for 
a given potential V{f). 

B. Partial Waves 

It is often useful to expand the scattering prob- 
lem in terms of partial waves, i.e. spherical harmonics 
Yi,mi {"^i f)- The result for the incoming plane wave is 

oo 

^ikz ^ i'j;(fcr)YioWV47r(2/ + l) , (3) 

where ji denotes the spherical Bessel function of order 
I. The outgoing wave is also written as a sum of partial 
waves. In this paper, we restrict the scattering problem 
to the case, where all outgoing partial waves have mi = 0. 
This is the case, e.g., if the potential is invariant under 
rotations aroimd the z-axis. The experiment in Ref. [16] 
does not have such a cylindrically symmetric potential, 
yet all outgoing partial waves have m; = 0, as we will see 
in Sec. VB. Hence, 

oo 

fid) - Yl mo{d)V^M2iTi) , (4) 

(=0 

where the partial- wave coefficients arc labeled With 
these expansions, one can reformulate the scattering 



problem. As a first step, one solves the scattering prob- 
lem for one incoming partial wave I' , where the asymp- 
totic form of the scattering state is 

-^f^Su'Yiom. (5) 

The outgoing partial waves have certain complex am- 
plitudes. These amplitudes form the so-called 5-matrix 
(or scattering matrix). Conservation of the number of 
particles implies that the S-matrix is unitary. In addi- 
tion, realistic Hamiltonians in atomic physics are invari- 
ant under time reversal, which implies that the 5-matrix 
is symmetric. 

The second step is to superpose the incoming partial 
waves with suitable amplitudes to obtain an incoming 
plane wave. This yields 

1 °° 1 21' -\- 1 

/' = ^EV27TT(^«'-'^"') ' 

where Sw is the Kronecker symbol. The task of calculat- 
ing /{{)) for a given potential F(r) is therefore replaced 
by the task of calculating the S'-matrix. Realistic poten- 
tials lead to selection rules for the angular momentum. 
Therefore, there are usually only few non- vanishing ma- 
trix elements in S. Calculating the ^-matrix is therefore 
often easier than calculating /(t?) directly. 

If more than one scattered partial wave is populated, 
then the differential cross section shows a spatial in- 
terference pattern between the partial waves. The to- 
tal cross section, however, shows no interference because 
/ r^foYjodri = 6w ■ Hence, 

OO OO 

1=0 1=0 

When calculating differential or total cross sections, spe- 
cial attention must be paid in the case of indistinguish- 
able particles, because then the two-particle wave func- 
tion needs proper symmetrization. We restrict the rest 
of this paper to the scattering of identical bosons, where 
the cross sections double for the even partial waves and 
vanish for the odd partial waves. 

C. Spherical Symmetry 

Things simplify if the potential is spherically symmet- 
ric. The quantum number I is then conserved, so that the 
S'-matrix is diagonal. Combined with unitarity this im- 
plies \Sii' I = 6ii'. All the information about the 5- matrix 
is therefore in the phases of the diagonal elements and 
one defines the scattering phase Si for each partial wave 
by 

Su = e^'^' . (8) 



3 



Note that the scattering phase is real and only defined 
modulo TT. The connection between /; and the ^-matrix 
simplifies to 



fi = 



Su-1 

2ik 



(9) 



For identical bosons, the total cross section for the l-th 
partial wave is 



(10) 



if I is even; and otherwise. Since 6i is real, ai has an 
upper bound 



^rnax = (2/ + 1) 



Sir 



'I -V 

which is called the unitarity limit. 



(11) 



the Schrodingcr equation must be solved along the radial 
coordinate. Various numerical methods have been devel- 
oped to solve this problem. If more than one collision 
channel is involved, then coupling between the channels 
must be taken into account. The corresponding calcula- 
tions arc called coupled-channels calculations. 

For atom-atom collisions (except for atomic hydrogen) 
there is another problem: ab-initio calculations for the 
interaction potentials arc not accurate enough to make 
realistic predictions for the cold-collision properties. As 
a solution, some quantities, such as the van-der-Waals 
coefficient Cf), arc treated as free fit parameters and ex- 
perimental input is used to constrain the model, in order 
to obtain realistic predictions for the cold-collision prop- 
erties (see e.g. Ref. [20]). 



III. RESONANCE SCATTERING 



D. Threshold Behavior 



S-Matrix 



The low-energy limit of the scattering properties is of- 
ten important. For reasons discussed in Sec. IV A, the 
energy where fc = is called the dissociation threshold. 
The behavior near threshold can often be expressed in 
terms of simple power laws. This was first systematically 
investigated by Wigner [17]. If the potential is spherically 
symmetric and has a long-range tail following a power 
law V{r) oc r~*, then one can show that (see p. 230 in 
Ref. [18]) 



fe^o r 0(fc2') if 2l<s- 
~ \ 0{k^-^) otherwise 



(12) 



Hence, Eqs. (7)-(9) imply that Si vanishes near threshold 

like k fi and that ai vanishes like f^. For s > 3 it follows 
that, in the low-energy limit, s-wave scattering dominates 

over all other partial waves and that Sq ~ 0{k). This 
motivates the definition of the s-wave scattering length 
[19] 



■ lim — 

/c->o k 



(13) 



The total scattering cross section for identical bosons is 
then 



's-»0 o 2 

a sna 



(14) 



The regime of cold collisions is characterized by ener- 
gies that are so low that only few partial waves have a 
noticeable scattering cross section. If only s-waves are 
important, the collisions are called ultracold. 



E. Coupled-Channels Calculations 

After performing the partial-wave expansion, the re- 
maining problem in scattering theory lies in the calcu- 
lation of the S'-matrix for a given potential. In essence. 



Resonance scattering relies on the presence of a quasi- 
bound state. A quasi-bound state is a discrete state just 
like a bound state, but with an energy above thresh- 
old. Hence, if one prepares population in this state, 
it will undergo spontaneous exponential decay into un- 
bound states. Examples for quasi-bound states are given 
in Sees. Ill E and IV. The decay of the quasi-bound state 
is a dissociation process, because a quasi-bound system 
decays into two unbound particles. The words decay and 
dissociation are therefore synonymous in the present con- 
text. 

In a scattering experiment, some fraction of the incom- 
ing flux can make the transition to the quasi-bound state 
and subsequently decay back into unbound states. Ob- 
viously, the probability to make this transition must de- 
pend on the energy difference between the incoming flux 
and the quasi-bound state. When the energies match, 
the population of the quasi-bound state is resonantly 
enhanced, while far-off resonance the population of the 
quasi-bound state becomes negligible. 

Near resonance, the S'-matrix is changed drastically. 
In this section, we consider only the case, where the S'- 
matrix is diagonal and only one partial wave has a reso- 
nance. With some effort one can show that the relevant 
S-matrix element is well approximated by a Breit- Wigner 
form [18, 21] (see appendix A for a derivation) 



Sii 



e''^' I 1 



E - Eres + ihT/2 



(15) 



Here Si ^ is the background value of the scattering phase 
for the relevant partial wave. This value is reached for 

energies far away from the resonance. Ej.es is the energy, 
at which the resonance occurs. The parameter F must be 
positive and can be interpreted as the decay rate of the 
quasi-bound state as discussed in Sec. IIIC. Generally, 
S'^^ and r can depend on E, but in the following we 
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partial wave, yielding 



with 



and 



bg {q + e-T" 



cot(5?» 



bg 



{2l + lf-^su.H\^ 



(19) 



(20) 



(21) 



Equation (19) is called a Beutler-Fano profile. This 
asymmetric resonance profile was first observed experi- 
mentally by Beutler in an autoionization experiment [22] 
and then explained by Fane [23, 24]. An example of a 
Beutler-Fano profile for g = 2 is shown in Fig. lb. 

A physical interpretation of the asymmetry in the 
Beutler-Fano profile is obtained easily when writing the 
^-matrix from Eq. (15) as the sum of a background scat- 
tered part S\? = e'^^^'i" and a resonantly scattered part 



FIG. 1: A typical scattering resonance, (a) Scattering phase 
calculated from Eqs. (17) and (18) for q — —cot5\^ = 2 
as a function of the dimensionless energy e. The horizontal 
lines show and + n, which are reached for e — »• ±oo. 
(b) Beutler-Fano profile for the cross section calculated from 
Eq. (19) for the same resonance. The horizontal line shows 
the background value. 



assume that they are independent of E within the width 
of the resonance. 

It is customary to introduce the dimensionless detun- 
ing of the energy from resonance 



Sii 



ore 
^11 



(22) 



e = ^(-E^ - Eres) ■ 



(16) 



Since \Sii \ = 1, one can again write Su = e^'*', where Si 
is real and one can easily show that 



with 



Si = 5f 9 + SI 



e = -cot5p* 



(17) 



(18) 



This means that 5i (E) has the form of an inverse tangent 
plus an offset, as shown in Fig. la. In particular. Si 
increases by tt, as the energy moves all the way through 
resonance. Right on resonance SJ^^ = 7r/2 and Si has an 
inflection point. 



B. Cross Section 

From the S'-matrix in Eq. (15), one can easily derive 
an expression for the total cross section of the relevant 



Note that Sfi'^" is usually not unitary. It is merely the 
resonant contribution to the S'-matrix, but not really an 
S'-matrix by itself. 

It is obvious from Eq. (15), that the phase of S'J";** 
changes by tt, when E moves all the way through res- 
onance. Hence the interference between S^^ and Sl{"^ 
changes from constructive to destructive, or vice versa. It 
follows from Eq. (19), that there is complete destructive 
interference at e = — q, while the unitarity limit Eq. (11) 
is reached at e = 1/q. At these points Si reaches and 
7r/2, respectively. The resonance position, i.e. e = 0, is at 
the inflection point of Si , which is usually not identical to 
the maximum of ct;. Note that q can be positive or neg- 
ative, so that the region of destructive interference can 
occur on either side of the resonance. For [^l S> 1 (i.e. 

~ 0), the Beutler-Fano profile is well approximated 
by a Lorentzian (except way out in the wings). 



C. Decay of the Quasi-Bound State 

Scattering theory is usually formulated as a time- 
independent problem. But in order to obtain a physical 
interpretation of the parameter F in Eq. (15), we now 
consider resonance scattering of a pulsed incoming wave. 
The incoming pulse can be expressed as a superposition 
of incoming plane waves, each of which can be described 
by time-independent scattering theory. We first consider 
the case where the energy width AE of the incoming 
pulse is large, i.e. AE ^ HT. In this case, one can show 
(see p. 254 in Ref. [18]) that the shape of the scattered 
wave packet has a tail at long times t. In the tail, the 
probability of detecting a scattered particle falls off like 
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FIG. 2: A quasi-bound state in the collision of two ^'^Rb 
atoms. Potentials for the s wave (dotted line) and the d 
wave (solid line) are shown as a function of radius. The cen- 
trifugal barrier for the d wave is clearly visible. In addition, 
the highest-lying s-wave bound state (horizontal dotted line) 
and the corresponding d-wavc state (horizontal solid line) are 
shown. The d-wave state lies above threshold (i.e. E > 0), 
so that it is only quasi-bound and gives rise to a scattering 
resonance. 

e~^*. This suggests the following interpretation: Some 
fraction of the population makes the transition to the 
quasi-bound state and decays from there with a rate T 
into outgoing waves. T is thus interpreted as the decay 
rate of the quasi-bound state. 

Since AE ^ hT, the minimum duration of the incom- 
ing pulse At ~ h/AE is much shorter than the mean 
lifetime of the quasi-bound state l/F and the tail can 
easily be distinguished from the background scattered 
wave packet. In the opposite limit, where AE <C hT, the 
S'-matrix is nearly constant within AE and the scatter- 
ing process leaves the temporal shape of the pulse almost 
unchanged. 

For time-independent elastic scattering, conservation 
of energy implies that the outgoing wave must have the 
same energy as the incoming wave. This is not so simple 
in resonance scattering of a short pulse. Here, the reso- 
nance plays the role of an energy filter. Only incoming 
energies with \E — Eres\ < ^T have a large probability to 
make the transition to the quasi-bound state. Only these 
energies are found in the exponential tail of the decaying 
wave function. 



D. Resonances near Threshold 



The energy dependence of F is particularly strong near 
threshold due to the centrifugal barrier. The centrifugal 
barrier arises from the centrifugal potential 

_ n'i{i + i) 

Vcentr- ^m.^dT^ ' ^''"^^ 

which depends on I and vanishes for the s wave. Compet- 
ing with the van-der-Waals potential — Cgr"^, the cen- 
trifugal potential dominates at large radius and the van- 
der-Waals potential dominates at shorter radius. This 
results in a centrifugal barrier of finite height, as shown 
in Fig. 2. Since the centrifugal potential vanishes at large 
radius, it does not shift the threshold energy. 

With the centrifugal barrier, the energy dependence 
of F near threshold can be understood as follows: For 
energies below the height of the centrifugal barrier, pop- 
ulation in a quasi-bound state must tunnel through the 
centrifugal barrier to decay into an outgoing wave. When 
the incoming energy approaches threshold, so does the 
outgoing energy because energy is conserved in time- 
independent elastic scattering. Thus when approaching 
threshold, the distance through which the particles must 
tunnel diverges. We thus expect that F vanishes near 
threshold, at least for 1^0. 

A more quantitative way to show that F must vanish 
near threshold is the following: When inserting the S'- 
matrix Eq. (15) into Eq. (9), one finds that /; cannot 
vanish faster than F/fc. Equation (12) for fi therefore 
sets a limit on how fast F must vanish. One can show 
that the following threshold law applies [25] 

F 0(fc2'+i) . (24) 

This is valid for all I, even if the potential has a long- 
range tail following a power law V{t) oc r^'^ [25]. Hence, 
for 2i > s — 3 we find that F vanishes even faster than 
kfi. For the decay into s- waves, the physical meaning of 
this threshold law arises from the density of final states, 
which is proportional to k (see e.g. Refs. [10, 26]). 

The Beutler-Fano profile for the cross section is typi- 
cally distorted for resonances near threshold by the en- 
ergy dependence of S'^^ and F, as well as by the factor 
in the unitarity limit Eq. (11). follows the usual 

threshold law S^^ '^^^ kfi with /; following Eq. (12). Note 
that for the molecule dissociation experiment in Sec. VI, 
k in Eq. (24) is the wave vector of the outgoing wave. 



While the consideration of pulsed scattering yields an 
intuitive interpretation of the parameter F in terms of 
a decay rate, one can show that for time-independent 
scattering the parameter F in Eq. (15) can depend on 
the collision energy. This is because, F is not merely a 
property of the quasi-bound state. Instead, it describes 
the decay into outgoing waves, so that the energy of the 
outgoing waves is important for F. We will now explain 
this in more detail. 



E. Shape Resonance in Rb 

A shape resonance is a scattering resonance, which is 
caused by a quasi-bound state behind some potential bar- 
rier. In this section we discuss a specific example, namely 
the d-wave shape resonance in the collision of two cold 
^^Rb atoms. We assume that both atoms are initially 
prepared in the lowest hyperfine state \f,mf) = |1,1) of 
the electronic ground state. This two-atom system has a 
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FIG. 3: Shape resonance for the scattering of two Rb atoms 
in the hyperfine state |1, 1) at S = 0. (a) Scattering phases 

and (b) cross sections arc shown for the s wave (dotted lines) 
and the d-wave (solid lines). The d-wave shape resonance is 
clearly visible. Other partial waves yield no noticeable scat- 
tering in this energy range. 



quasi-bound d-wave state localized behind the centrifu- 
gal barrier as shown in Fig. 2. Also shown in Fig. 2 are 
the s-wavc and rf-wavc potentials and the highest-lying 
s-wave bound state. The centrifugal potential raises the 
energy of the corresponding rf-wave state. In ^''Rb the 
energy of this state lies below the top of the centrifu- 
gal barrier, so that it is not pushed out into the uncon- 
strained continuum. But at the same time it lies above 
threshold and can therefore decay into the continuum 
by tunneling through the centrifugal barrier, so that it 
is a quasi-bound state. Tunnc;ling also allows incoming 
flux in a scattering experiment to populate the state. 
Since the terms in the interaction Hamiltonian that can 
change / arc weak (sec Sec. IV C), this quasi-bound d- 
wave state couples almost exclusively to the d-wave con- 
tinuum, so that a shape resonance is caused for incoming 
d waves, but no noticeable resonance is caused for incom- 
ing s waves. 

The partial-wave scattering phases and cross sec- 
tions determined from a coupled-channels calculation are 
shown in Fig. 3 for zero magnetic field, B = 0. The 
threshold law Eq. (12) for the s wave predicts that 

Sq ^^'^ 0(i?^/^) and that ctq approaches some finite value 
for small E. This is the case for small values of E, as is 
clearly seen in Fig. 3. 

The d-wavc scattering phase and the d-wavc cross sec- 
tion both show clear signatures of the shape resonance. 



In contrast to the general discussion above, decreases 
considerably over the width of the resonance. As a con- 
sequence, 62 increases by noticeably less than tt. Prom 
the threshold law Eq. (12) with s = 6, one expects 

^~ 0{E^), which leads to a strong asymmetry of the 
resonance in the cross section merely from the change in 
the background value. The additional inherent asymme- 
try of the Beutler-Fano profile is small for this resonance 
because \q\ ^ 1. 

The position and width of the shape resonance can be 
read-off from Fig. 3a. The inflection point of 62 (E) lies 
at Eshape = fc_B X 273 (jK. According to Eq. (18), the 
slope at the inflection point is d52/dE = 2/(?i,r), which 
yields F = 17 MHz. 

While a coupled-channels calculation is needed to ob- 
tain Fig. 3, one can already make a rough estimate for F 
if only Eshape and Cq = 4707 a.u. [27] are known (1 a.u. 
= 9.573 X 10"^" Jm*'). To this end, we approximate the 
potential as = Vcentr ^ C^r^^, which is a good ap- 
proximation at large radius. With this, we first calculate 
the tunneling probability in the WKB approximation for 
a spherical wave at energy Eshape and obtain 25%. Al- 
ternatively, the tunneling probability can be calculated 
analytically using a near-threshold approximation as dis- 
cussed in Ref. [28]. Equation (57) in Ref. [28] yields 23%. 
Second, we consider a classical particle with mass m^ed 
released at rest at the outer classical turning point of the 
quasi-bound state. In the potential V = Vcentr — Cer~^, 
this particle will roll down to r = in 6.3 ns. The ap- 
proximative treatment of the potential is reasonable, be- 
cause most of the time is spent near the outer turning 
point, where V = Vcentr — Cer~^ is a good approxima- 
tion. We assume that the particle is simply reflected at 
r = so that the round-trip time is 12.6 ns. We thus ob- 
tain F ~ 0.25/(12.6 ns) = 20 MHz, which is quite close 
to the above result. 

All other hyperfine levels of the electronic ground state 
of ^''Rb have a similar shape resonance with almost the 
same energy and lifetime. The first experimental ob- 
servation of this quasi-bound state had been made in a 
photoassociation experiment [29] . For the hyperflne state 
]2,2), this shape resonance was recently investigated in 
two scattering experiments [30, 31]. Other atomic species 
also have shape resonances in cold collisions, as measured 
e.g. in Refs. [32-38]. 



IV. MAGNETICALLY TUNABLE FESHBACH 
RESONANCES 

A. General 

In many scattering experiments, the particles have in- 
ternal degrees of freedom, such as spin. For each spin 
state, there is a different potential. These potentials are 
referred to as the scattering channels. We consider the 
situation sketched in Fig. 4, where we picked two chan- 
nels with different threshold energies. We assume that 
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FIG. 4: Scheme of a Feshbach resonance. The interaction 
Hamiltonian can cause transitions between the closed channel 
and the incoming flux in the open channel. A Feshbach res- 
onance occurs when the energy of the incoming flux matches 
the energy of a closed-channel quasi-bound state. 



the incoming flux has the spin configuration of the lower 
channel and that the energy of the incoming flux is below 
the threshold of the upper channel. In this situation, no 
flux can emerge in the upper channel for energetic rea- 
sons. This channel is therefore energetically closed, while 
the entrance channel is always open. If the entrance chan- 
nel is the only open channel, then a two-body collision 
cannot be inelastic, i.e. the spin states before and after 
scattering must be identical. 

A Feshbach resonance arises if incoming flux in the 
open channel is resonant with the energy of a bound state 
in a closed channel. We call this bound state involved 
in the Feshbach resonance the "molecular state". For 
the resonance to occur, the interaction Hamiltonian must 
be able to flip the spins, in order to cause transitions 
between the two channels. These spin flips lead to decay 
of the moleciilar state into unbound open-channel states, 
so that the molecular state is only quasi-bound. 

The potentials shown in Fig. 4 schematically represent 

Born-Oppcnheimcr potentials for the collision of two al- 
kali atoms. The potentials have a long-range van-der- 
Waals tail, a deeply bound region due to the exchange 
interaction, and a repulsive part at very short radius due 
to the Coulomb interaction of the nuclei and a repulsive 
exchange interaction of the overlapping electron clouds. 
The energies in Fig. 4 are not to scale. For ^^Rb, for ex- 
ample, the dissociation energy of the singlet potential is 
k-B X 5750 K [39], whereas the hypcrfinc splitting between 
the different thresholds is ks x 0.33 K. Each potential 
has many bound states. For clarity, only one of them is 
shown in Fig. 4. 

The general treatment of scattering resonances given in 
Sec. Ill is also applicable to Feshbach resonances. In par- 
ticular, one can show that the S'-matrix is again given by 
Eq. (15). To obtain this result, one often uses a formal- 
ism based on Green's functions and projection operators 
introduced by Feshbach [40, 41], in order to clearly sepa- 
rate the open-channel and closed-channel subspaces. For 



a discussion of this formalism in the context of ultracold 
gases see e.g. Ref. [26]. 

In the collisions of cold alkali atoms, there are typ- 
ically many Feshbach resonances. The molecular state 
is a vibrationally highly excited state. For homonuclear 
inok;cules, the radiative decay rates into lower vibrational 
states are negligible due to lack of an electric dipole mo- 
ment. But in a dense sample, inelastic collisions of the 
molecule with another atom or molecule can lead to sig- 
nificant rates for vibrational de-excitation. These pro- 
cesses are neglected in this paper. 

Feshbach resonances in the collisions of cold atoms dif- 
fer substantially from most scattering resonances in other 
fields of physics, insofar as the energy of the molecular 
state can be tuned noticeably by applying an external 
magnetic field B. The tunability arises from the differ- 
ent magnetic moments of the different channels. If this 
difference is, say, one Bohr magneton, a magnetic field 
of 1000 G creates a shift of ks x 0.07 K. Therefore low 
temperatures of the incoming flux arc needed, otherwise 
thermal broadening would render this shift insigniflcant. 

Let Bres denote the magnetic field, at which the energy 
of the moleciilar state matches the open-channel thresh- 
old. Then the magnetic-field dependence of the energy 
Eres, at which the Feshbach resonance occurs, can be 
approximated linearly by 

Eres{B) = {B - Bres)^l^ (25) 

as long as [S — -Bres I is not too large. denotes the 
difference in the magnetic moments of the two channels. 

can be positive or negative. If B approaches Bres, 
then Eres ^ i.e. the Feshbach resonance occurs at en- 
ergy zero. A small correction to Eq. (25) is necessary if 
B is very close to Bres ■ This correction is typically of the 
order of AB (see below). This so-called resonance shift 
is neglected in the discussion in the present paper. But 
it is included in our coupled- channels calculations. 

For Eres{B) < the molecular state is truly bound, 
which means that the state cannot decay. When varying 
B so that Eres increases, dissociation abruptly sets in 
at Eres{B) = 0. Therefore the low-energy edge of the 
continuum is called dissociation threshold. 

An interesting situation arises, if one ramps B through 
Bres in the direction such that Eres moves from above 
to below threshold. This converts the molecular state 
from quasi-bound to truly bound. If incoming atomic 
flux is present during the ramp, then population that was 
transiently in the molecular state during scattering while 
above threshold, will remain in the state after the thresh- 
old is crossed. Thus long-lived molecules can be pro- 
duced. These can later be dissociated at will by ramping 
B back through Bres- This technique was used in several 
recent experiments to produce ultracold molecules from 
ultracold atomic gases [1 9]. Between the association and 
the dissociation, the atomic and the molecular cloud can 
be spatially separated by applying a Stern-Gerlach field. 

We emphasize the difference between a Feshbach res- 
onance and a shape resonance. A shape resonance oc- 
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curs in single-channel scattering. The quasi-bound state 
is typically localized behind the centrifugal barrier. It 
could theoretically become a truly bound state by in- 
creasing the potential depth, but experimentally, one can 
usually not tune its energy. A Feshbach resonance is a 
multi-channel resonance. The quasi-bound state has a 
spin configuration other than the incoming flux. The 
energy of the quasi-bound state can be tuned with a 
magnetic field and this state can become a truly bound 
state by tuning its energy below the open-channel thresh- 
old. Despite their differences, both types of resonances 
have in common that there is a quasi-bound state above 
threshold. Moreover, both types of resonances are well 
described by the same Breit-Wigner expression for the 
^-matrix Eq. (15). 



B. Low-Energy Feshbach Resonances 

As discussed in Sec. II D, s-wave scattering usually 
dominates at low energies. If the magnetic field is held 
near a Feshbach resonance at Bres, then the low-energy 
scattering is affected, of course. According to Eqs. (12) 
and (24), 6o and F both vanish like 0{k) for A: ^ 0. 
One can thus linearize the tangent in Eq. (18), yielding 
(5(7" - tan = nr/[2{Eres-E)]. The definition of the 
scattering length Eq. (13) then yields 



abg 



^ lim — . 



2Eres *:->0 k 



(26) 



Since F ~" 0{k), this expression is well defined. Insert- 
ing Eq. (25) for Eres one obtains 



a = ttftg I 1 - 



AB 



B — Br 



(27) 



where the (magnetic-field) width of the Feshbach reso- 
nance is defined as 



lim — . 



2abgA/i fe^o k 



(28) 



F is always positive, while A/x and a^g can be positive or 
negative, independent of each other. Therefore AB can 
be positive or negative. 

Low-energy scattering in the vicinity of a Feshbach res- 
onance can thus be described analytically with only three 
parameters a^g, Bres, AS. The resulting pole and zero in 
the total cross section a = 8na^ result from interference 
between background scattered and resonantly scattered 
wave, just like in the Beutler-Fano profile. Various ex- 
periments [42-48] with ultracold atoms observed the be- 
havior predicted by Eq. (27). 

For sufficiently small k, Eq. (28) can be solved for F. 
Inserting E = /i^fc^/(2mred), one obtains the threshold 
law for the decay rate 



^ 2ABA/X /- 

1 = — ^2 abgy/2mredE 



(29) 



Recently, this behavior was also experimentally observed 
with ultracold atoms [10, 11]. 



C. Selection Rules 

As mentioned above, one requirement for a Feshbach 
resonance is that the interaction Hamiltonian must be 
able to flip the spins to make transitions between the 
two relevant channels. This section deals with this issue 
in more detail. 

At large radius, the spins of two colliding ground-state 
alkali atoms are specifled in terms of the hyperflne quan- 
tum numbers |/i,m/i) and |/2,?ti/2) of the two atoms. 
Together with Z, to;, one obtains a complete set of quan- 
tum numbers. The atomic hyperflne spins can be added, 
yielding the total spin F = f\ + f^- The corresponding 
quantum numbers are F, mp ■ 

At shorter radius, the exchange interaction Vex is the 
dominant term in the interaction Hamiltonian, so that 
the spins of the valence electrons are coupled to a to- 
tal electronic spin 5"; and the singlet {S = 0) and triplet 
(5* = 1) potentials differ drastically. Hence, the hyperflne 
quantum numbers /i, m/i, /2, my^2 are not good quantum 
numbers at short radius. When writing Vpx as a matrix 
in the hyperflne basis, it therefore has large off-diagonal 
elements, which means that transitions between differ- 
ent hyperflne states are possible. Vex is spherically sym- 
metric and thus conserves Lto,;. For incoming s waves. 
Vex can therefore cause Feshbach resonances only if the 
molecular state is an s-wave state. Since Vex creates only 
forces internal to the system, the total angular momen- 
tum / -|- _F is conserved. Since to; is conserved, mp is 
conserved, too. 

In addition, there are much weaker terms in the in- 
teraction Hamiltonian. The strongest of these terms is 
the spin-spin interaction Vss, which is the sum of the 
magnetic dipole-dipole interaction of the valence elec- 
trons and the second-order spin-orbit interaction for the 
valence electrons. Vis can change Z,to; because it is 
not invariant under spatial rotations. It causes transi- 
tions according to the selection rules Ai = or ±2 and 
|Ato;| < 2. For incoming s waves, Vss can therefore 
cause Feshbach resonances for d-wave molecular states. 
Since Vss is much weaker than Vex, the resulting inter- 
channel coupling is typically also much weaker. Hence, 
Feshbach resonances caused by Vss are usually much nar- 
rower than those caused by Vex- Since Vss creates only 
internal forces, the total angular momentum l+F is again 
conserved. 

Even weaker terms in the Hamiltonian can cause other 
narrow Feshbach resonances, such as in ^^^Cs near 20 
Gauss, where an incoming s wave is coupled to a g-wave 
molecular state [49]. A^ is always even for the Feshbach 
resonances in atomic collisions, because the interaction 
Hamiltonian conserves parity. The only fundamental in- 
teraction that does not conserve parity is the weak inter- 
action, but that is negligible here. 
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When an external magnetic field B is applied, the total 
angular momentum l + F is no longer conserved, because 
the external field creates external forces. We consider 
only the case where B points along the z axis, so that 
rotational symmetry around the z axis implies the con- 
servation of nil + TTiF- Note that if the magnetic field 
is strong, /i, /2 are no longer good quantum numbers at 
large radius. 



V. RESONANCES WITH MANY PARTIAL 
WAVES 

A. S-Matrix 

In this section, we consider the case where a quasi- 
bound state couples to unbound states in more than one 
partial wave. Again, a Breit-Wigner form is obtained for 
the S-matrix (see p. 411 in Ref. [18]) 



s = s^ni 



iA 



E - Eres + ihT/2 



(30) 



where S''^ and A are matrices and 1 is the identity ma- 
trix. While this is often discussed in the context of res- 
onances with couplings between different channels, we 
here use it for resonances that involve different par- 
tial waves. We assume that S^^ is diagonal and define 

Til, = e'^^'''~^''i' '>Aii,/h. We thus obtain 



Sw 



ihTw 



E — Er 



ihT/2 



(31) 



Note that this reduces to Eq. (15) in the single partial- 
wave case, where Fn' = TSu'Suo and denotes the one 
partial wave that has a resonance. 

We assume that there is only one quasi-bound state 
that causes the resonance. Hence, the matrix A is of 
rank 1 (sec p. 406 in Ref. [18]). Combined with unitarity 
and symmetry of the 5'-matrix, this implies that all F;;/ 
are real and that they fulfill 



(32) 
(33) 



We will see in Eq. (43) that the decay rate into the l- 
th partial wave is given by F;;. This quantity cannot 
be negative. The threshold law Eq. (24) applies to each 
diagonal element F;(. The total decay rate F is the sum 
of the partial decay rates. 

This situation explicitly allows coupling between dif- 
ferent partial waves, so that \Su'\ ^ Sw. In particular, 
\Sii\ ^ 1 and the scattering phases 5i can be complex. 
The S'-matrix is still unitary and the number of parti- 
cles is conserved. But the unitarity limit Eq. (11) for ai 
known from spherically symmetric potentials can be ex- 
ceeded here, because flux can be redistributed between 
partial waves. 



When calculating the partial-wave components of the 
total cross section, one finds an expression that closely 
resembles a Beutler-Fano profile 

^ ^bg (e + Re{Q})^^+(l+Im{Q})^ ^^^^ 

with e and a^^ from Eqs. (16) and (21) and with the 
complex number 



Q 



1 V- ^^sy Tll'VW+l 



(35) 



While Q is a lengthy expression, the key result is that 

it is independent of e. Hence, cr; is a fairly simple func- 
tion of energy, namely the product of a parabola and 
a Lorentzian, just like for the normal Beutler-Fano pro- 
file. The full width at half maximum (FWHM) of the 
Lorentzian as a function of E is hT just like in the single 
partial-wave case. But the minimum of ai can be above 
zero and the maximum can be above or below the unitar- 
ity limit. In the single partial- wave case {Tw = rSu'Sug), 
one obtains Q = {—i — cot 5'l^)5iig and Eq. (34) reduces 
to the normal Beutler-Fano profile Eq. (19). 



B. Combination of a Shape Resonance and a 
Feshbach Resonance in ^^Rb 

We apply this formalism to ®^Rb with both atoms en- 
tering in the hyperfine state |1, 1). This is the absolute 
ground state for atomic ^^Rb so that inelastic two-body 
collisions cannot occur. This system has a Feshbach res- 
onance at Bres = 632.45 G [27] with AB = 1.3 mG 
[11]. The corresponding molecular state is a d-wave 
state, which is coupled to incoming s, d, and g waves 
by the spin-spin interaction (see Sec. IV C). For the en- 
ergy range considered here, g-wa,ve scattering is negligi- 
ble. Weaker terms in the Hamiltonian can couple to even 
higher partial waves, but that is also negligible here. 

Figure 5 shows the s and d-wave components of the 
total cross section, tJo and CT2 with mi = (see below), 
for different values of the magnetic field B. These results 
were obtained from a coupled-channels calculation. One 
can clearly see the narrow Feshbach resonance sitting on 
top of the background, which is modulated due to the 
broad shape resonance already known from Fig. 3. The 
energy Eres, at which the Feshbach resonance occurs, can 
be tuned with the magnetic field according to Eq. (25). 

F in Eq. (31) denotes the total decay rate of the molec- 
ular state. The decay rate of the other quasi-bound 
state that causes the shape resonance does not explic- 
itly occur in Eq. (31), because the open-channel physics 
including the shape resonance is contained in the energy- 
dependence of S'^^ and F;;/ . 

When moving Eres{B) through the shape resonance, 
the form of the cross sections ai near the Feshbach reso- 
nance changes, as seen in Fig. 5. For Eres{B) < Eghape, 
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FIG. 5: Combination of a shape resonance and a Feslibach 
resonance for scattering of *^Rb in state jl, 1). Tlic partial- 
wave components cr; of tlie total cross section are shown for 
the s wave (dotted line) and the d wave (solid line). The 
magnetic field B is held at various values above the Feshbach 
resonance at Bres ~ 632 G. The Feshbach resonance is much 
narrower than the shape resonance. By changing B, the po- 
sition of the Feshbach resonance can be tuned through the 
shape resonance. 



the Feshbach resonance increases the cross sections on 
the low-energy side of the Feshbach resonance and re- 
duces the cross sections on the high-energy side. For 
Eres{B) > Eshape this is reversed. This is because when 
moving E through the shape resonance, 82' changes by 
almost TT, as seen in Fig. 3a. This phenomenon is called 
g-reversal, see e.g. Ref. [50]. 

We now discuss, why m; = for all outgoing partial 
waves. As mentioned in Sec. IV C, rotational symme- 
try around B, which points along the z axis, implies 
that mi + mp is conserved. Since the incoming chan- 
nel has rrii;' = 2, m; = 0, all outgoing channels must have 
mp + m,i =2. For energetic reasons, mp = 2 is the only 
possible spin state for the outgoing flux, so that all outgo- 
ing waves must have mp = 2, mi = 0. Note that during 
the collision, couplings to all states with mp + mi = 2 are 
possible. This is a necessary ingredient, since the molec- 
ular state is an almost pure mp = A,mi = —2 state. 

Note that in the case of the above-mentioned reso- 
nance near 632 G, Vgs is required for the molecules to 



decay, since to/ has to change by +2 for both the out- 
going s and d wave, leading to comparable amplitudes 
for both partial waves. If the molecular state were an 
s-wave state instead, Vss would still couple it to the out- 
going d wave, but the much stronger Vex would create 
a strong coupling to the outgoing s wave, resulting in a 
very small outgoing d-wave fraction. Conversely, if the 
molecules had I = 2,mi = they could decay into the 
outgoing d wave by the strong Vex since I and to; would 
not need to change. For decay into the outgoing s wave, 
the much weaker Vgs would be needed as I would need 
to change by —2 and the population of this partial wave 
would therefore be strongly suppressed as compared to 
population in the li-wave. Such resonances do exist, e.g. 
in s^Rb at 551.47 G and 831.29 G [27]. Unfortunately, 
these two resonances are so narrow {AB ~ 0.2 mG each), 
that creating molecules at these resonances is difficult. 



C. Extracting the Partial Decay Rates 

The 5-matrix for the above situation was numerically 
calculated on a fine grid in the i?-i3-plane, in order to 
extract the decay-rate matrix F;;/ and the background 
scattering-phases 6^^ as defined in Eq. (31). The extrac- 
tion of the ^J"^ is easy, because they contain only the 
open-channel physics. This does include the shape reso- 
nance, so that the S'l^ depend on energy, but it does not 
include the Feshbach resonance, so that the 6^^ depend 
hardly on the magnetic field. The values of S'l^{E) can 
therefore simply be read off from the 5-matrix for pretty 
much any B sufficiently far away from the Feshbach res- 
onance. Since the open-channel physics is almost inde- 
pendent of B, the result for d'^^ is essentially the same as 
shown for i? = in Fig. 3a. 

As discussed in the context of Eq. (24), the decay-rate 
matrix of the molecular state Fo' depends on energy. But 
it depends hardly on the magnetic field. Hence, the ex- 
traction of F;;' (E) from the numerical results for S{E, B) 
is also fairly easy. To this end, we take the modulus 
squared of Eq. (31) for Z 7^ V and insert Eq.(25), yielding 



\Siv 



12 li^l' 



n'Tl,{E)/At,^ 



{-B + Bres + E/Anf + h^T^{E)/{2An)^ 



.(36) 



Considered as a function of B at constant E, this is sim- 
ply a Lorentzian. We fit this to the numerical results 
for \So2{E, B)]"^ at constant E and obtain three fit pa- 
rameters |Fo2(-E)/A^|, |F(£;)/A^|, and {Bres +E/Aii). 
A combination of the results of the last fit parame- 
ter for various values of E yields Bres = 632.3 G and 
A/i = kB X- 224 fiK/G = 3.33 with the Bohr mag- 
neton j2B- The deviation between the theoretical and 
experimental value for Bres is no problem, as long as the 
comparison between experiment and theory is performed 
in terms oi B — Bres ■ 
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FIG. 6: Decay rates of the molecular state for the 632-G 
Feshbach resonance in ^'^Rb. The partial decay rates into the 
s wave (dotted line) and d wave (solid line) are shown. The 
fi-wavc shape resonance obviously has a drastic effect on the 
d-wave decay rate. 

Knowing A^, one obtains |ro2(i?)| and T{E). Using 
Eqs. (32) and (33), one easily obtains Toq{E) and T22{E), 
except for the ambiguity of which is which. This can 
easily be resolved by inspection of the diagonal elements 
of S{E,B). In addition, inspection of the off-diagonal 
elements of S{E, B) yields the sign of ro2(-E'), which is 

= = —\/TooT22 (37) 

for all E in the present calculation. Thus, the complete 

decay-rate matrix F^' {E) is extracted. The partial decay 
rates Too {E) and r22 {E) of the molecular state are shown 
in Fig. 6. 

In order to demonstrate the quality of the fit, the re- 
sults of the coupled-channels calculation are compared to 
the fit curves in Fig. 7. Parts (a) and (b) show modulus 
squared and phase of the S'-matrix elements, respectively. 
For clarity, some quantities that are trivially related to 
the others arc not shown. Symmetry and unitarity of the 
S matrix imply ^20 = 5*02, |5'ooP = \S22? = 1 - \So2? , 
and argS'02 = (7r-|-arg6'oo + arg6'22)/2. All curves shown 
in Fig. 7 are well described by only four parameters 5^^ , 
^'2 1 ^02, and r. The excellent agreement between the fit 
and the coupled-channels results demonstrates that the 
Breit-Wigner form Eq. (31) is a very good approxima- 
tion. 



VI. MOLECULE DISSOCIATION 
A. General 

In this section, we discuss the decay of the molecular 
state as observed in Ref. [16]. As already described at 
the end of Sec. IV A, molecules are formed by ramping 
B through B^es ■ A Stern-Gerlach field then removes the 

incoming flux from the spatial region of interest and fi- 
nally the molecules are dissociated by jumping B above 




633.35 633.40 633.45 633.50 633.55 
e(G) 

FIG. 7: Fitting to the 5-rriatrix elements. Results from the 
coupled-channels calculation (circles) are shown versus mag- 
netic field B for a fixed energy of fcs x 255 /uK. The fit curves 
(solid lines) are hardly visible, because they agree so well with 
the coupled-channels results. 



Br-es and holding it there at a fixed value. For further 
experimental details, see Refs. [3, 16]. 

Due to conservation of energy during the decay, the 
mean energy of the outgoing wave is given by 

E = Eres{B) (38) 

with Eres from Eq. (25). Due to the finite lifetime of the 
molecular state, the energy of the outgoing wave has a 
width of ^r. 

Note the difference between this dissociation experi- 
ment and time-independent elastic scattering: for time- 
independent scattering, the energy of the outgoing wave 
must be identical to the energy of the incoming wave, 
whereas in the decay experiment described here the en- 
ergy of the outgoing wave is adjusted with B. In the 
decay experiment, the energy of the outgoing wave of up 
to fcs X 500 /iK is typically much larger than the en- 
ergy of the incoming wave with E < ks ^ ^ ^K, from 
which the molecules were originally formed. This is pos- 
sible because the magnetic-field ramp makes the Hamil- 
tonian explicitly time dependent, so that energy is not 
conserved in a simple-minded fashion. Ultimately, the 
additional energy in the outgoing wave comes from the 
power supply that creates the magnetic-field ramp. 

In the decay experiment, the original incoming wave 
has such low energy, that only incoming s-waves are rel- 
evant. This implies that m; = for the incoming wave. 
As discussed in Sec. VB, this implies that all outgoing 
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partial waves must also have m; = 0. And again, only 
s and d waves are important for the dissociation in the 
energy range considered here. 



B. Connection between Scattering and Decay 

The decay experiment is closely related to the scat- 
tering experiment discussed in Sees. VB and VC. The 
link can be established by inspection of the scattering 
wave function, which is obtained from the combination 
of Eqs. (5) and (22) 



i-iy 



-ikr 



ikr 



ikr °° 



-V E^«'^'oW- — E^«''^'oW . (39) 



/=0 



1=0 



The first term is the incoming wave, the second term 
is the background scattered wave and the third term is 
the resonantly scattered wave. The resonantly scattered 
part consists of population that made the transition to 
the molecular state and then decayed back to the open 
channel. In the decay experiment in Ref. [16], the Stern- 
Gcrlach field removed the incoming wave and along with 
it the background scattered wave. Hence, these two terms 
must be removed to describe the decay experiment 



Jkr °° 



(40) 



1=0 



This crucial step makes the connection between scatter- 
ing and dissociation. 

Prom Eq. (31) one obtains 



E - Eres + 



(41) 



Inserting this and Eq. (37) into Eq. (40) and choosing 
/' = yields 



We abbreviate 



oikr 



E - Eres + 



/f^y2o(^)) • (42) 



PI - 

5rel = (^2^ - 



(43) 
(44) 



and obtain 
i^decay{r) g{r, E)(^^/p^Yoo - e 



i6r 



fi2Y2o{'d)) (45) 



where g{r,E) is the radial wave function. The time- 
domain version thereof g{r, t) is related to g{r, E) by a 
Fourier transform. 



Equation (33) implies that /3o + /?2 = 1 so that /?; is the 
branching ratio for decay into the ^-th partial wave. Srei 
is the relative phase between the two partial waves. The 
above definitions of /?; and Srei do not explicitly depend 
on E. But they do involve TwiE) and S'^^E). The 
energy-dependence of these quantities within the width 
of the resonance hV is negligible, so that one can simply 
evaluate these quantities at E = Eres{B). 

In conclusion, we showed how the dissociation of ultra- 
cold molecules into more than one partial wave is related 
to a scattering experiment. This makes it possible to use 
coupled-channels calculations for scattering experiments 
to analyze dissociation experiments, such as the one in 
Ref. [16]. The analysis in Ref. [16] shows that the exper- 
iment agrees well with the theory described here. 
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APPENDIX A: POLES OF THE S-MATRIX 
1. Analytic Continuation 

Insight into scattering resonances can be gained in a 
very general formalism that makes no use of the specific 
form of the potential. The starting point is the observa- 
tion that the S-matrix is typically an analytic function of 
k. The key idea is then to consider the analytic contin- 
uation of S{k) into the complex fc-plane. One can show 
that this continuation is unique, but cannot always cover 
the whole complex plane. The physical meaning of k as 
the magnitude of the wave vector requires it to be real 
and non-negative. But the continuation into the com- 
plex plane will offer additional physical insight, as we 
will see in the following. For simplicity, we consider only 
a spherically symmetric single-channel potential in this 
appendix. More details about the topics discussed in this 
appendix can be found in Ref. [18]. 

In this discussion, one usually considers the Jost func- 
tion fi{k) instead of the ^-matrix. Like the S'-matrix, 
the Jost function is also defined by the coefficients in the 
scattering state Eq. (5), but with a different normaliza- 
tion 

(p—ikr pikr \ 

i-iyMk)-^-fnk*)—)YioW (Al) 

The Jost function also has a unique analytic continuation 
into the complex k plane. For most potentials, fi{k) is 

analytic everywhere in this plane, except for the negative 
imaginary axis. Note that fi{k) and fi{k*) cannot both 
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vanish for the same value of k. As a consequence, fi{k) 
cannot vanish if k is real. 

The S-matrix is easily obtained from the Jost function: 



Sw{k) 



~m 



6u' . 



(A2) 



The Kronecker symbol comes from our assumption of 



spherical symmetry. If k is real, this implies \Sii 



i.e. unitarity. 

Note that the S'-matrix has a pole at k if fi{k) = 0. 
In other words, poles of the ^-matrix correspond to ze- 
ros of the Jost function. If k is real, the Jost function 
cannot vanish, so that all zeros of the Jost function must 
lie either in the upper half-plane, i.e. Im{fc} > 0, or in 
the lower half-plane, i.e. Im{fc} < 0. We will now discuss 
what the physical meaning of the zeros of the Jost func- 
tion is. This depends on the half plane, in which the zero 
is. 



only those poles that lie close to the positive real axis and 
when ignoring some rather special cases (see p. 241 in 
Ref. [18]), one can think of this correspondence as being 
one-to-one. 

The linear approximation further below is usually 
made in terms of energy, rather than k. Obviously, one 
can substitute E = Ti? k"^ / {2mred) and obtain fi{E). Care 
must be taken, because the transition from A; to £^ is a 
two-to-one mapping, so that fi{E) is a function on a two- 
sheeted Riemann surface. Im{A;} > corresponds to the 
first sheet of E (also called physical sheet) and Im{A;} < 
corresponds to the second (or unphysical) sheet of E. 

If fi{ko) = with ko in the lower half-plane, then fi{E) 
has a zero on the unphysical sheet at Eq = ?t^kQ/{2mred)- 
We assume that fi{E) has a simple zero at Eq, so that 
near Eq we can use a linear approximation 



fi{E) ^ c{E - Eo) 



(A3) 



2. Bound States 

For complex values of fc, Eq. (Al) has one term that 
increases exponentially for large r, and one term that 
decreases exponentially. We assume that there is a 
point fco ill the upper half-plane with fi{ka) = 0. Here, 
Eq. (Al) has only one term. This term falls off expo- 
nentially, because we assumed that fco is in the upper 
half-plane. Hence. Eq. (Al) is a normalizable eigenstate 
of the Hamiltonian. In other words, this is a bound state. 

The bound-state energy is h"^ k^ / {2mred) ■ Since the 
Hamiltonian is hermitian, this energy must be real. Since 
fco is in the upper half-plane, it follows that fco ~ ia 
with a real and positive. The bound-state energy is then 
—h^a^/{2mred)- Conversely, one can show that if the 
Hamiltonian has a bound state with angular momentum 
/ and energy —ti^a'^/{2mred), then fi{ia) = 0. 

To summarize, bound states have a one-to-one corre- 
spondence to poles of the S-matrix in the upper half of 
the complex fc plane. All these poles must lie on the 
imaginary axis. 



3. Resonances 

We now investigate what happens if the Jost function 
vanishes at a point fco in the lower half-plane. Again, the 
first term in Eq. (Al) vanishes, but this docs not result 
in a normalizable eigenstate, because the remaining term 
increases exponentially for large r. Unlike before, fco does 
not have to lie on the imaginary axis. 

The Jost function can thus have a large number of 
zeros everywhere in the lower half-plane. Most of them 
are usually uninteresting. The only interesting ones are 
those, which lie close to the positive real axis because 
they can create resonances, as we will show now. The 
correspondence between resonances and zeros of the Jost 
function in the lower half-plane is not so clear-cut one-to- 
one as in the case of bound states. Still, when considering 



with a nonzero value of c = ^j^.^- We already men- 
tioned that fco must lie just slightly below the positive 
real fc axis to create a noticeable resonance. This implies 
that Eq also lies just slightly below the positive real E 
axis. Hence, there is some interval on the positive real 
E axis, where Eq. (A3) is a good approximation. From 
Eq. (A2), one obtains 



Sii 



c*{E-Eo)* 
c{E-Eo) 



(A4) 



for E in this interval on the real axis. We abbreviate the 
prefactor that is independent of E as 



and we split Eq into its real and imaginary part 



(A5) 



(A6) 



Here, Eres and F must be positive, because we assumed 
that the zero of the Jost function lies slightly below the 
positive real fc axis. We thus obtain the Breit-Wigner 
expression 



ihv 



E - Eres+ifl^/2 



(A7) 



The real and imaginary parts of Eq are thus identified as 
the position and the width of the resonance. 

With this approach to the Broit-Wignor form, Sff^ and 
F are independent of E. If the resonance is narrow, i.e. 
F is small, then this is a good approximation. Broader 
resonances can be included in the formalism by allowing 
S^^ and F to depend on E. For a very broad resonance, 
the background scattering phase can change considerably 
over the width of the resonance, and it becomes question- 
able whether one really should regard this as a resonance. 
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